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Abstract. Current quantum cascade lasers based upon conduction band 
electron transitions are predominantly TM (electrical field normal to the epitaxial 
direction) polarized. Here we present a study of localized defect modes, with 
the requisite TM polarization, in connected square and hexagonal lattice two- 
dimensional (2D) photonic crystals for application as quantum cascade laser 
resonators. A simple group-theory based analysis is used to produce an 
approximate description of the resonant modes supported by defects situated at 
different high symmetry points within the host photonic lattice. The results 
of this analysis are compared with 2D finite-difference time-domain (FDTD) 
simulations, showing a close correspondence between the two analyses, and 
potential applications of the analysis in quantum cascade laser design are 
considered. 



PACS numbers: 42.70.Qs, 42.55.Sa, 42.60.Da, 42.55.Px 

1. Introduction 

Optical microcavities have been the subject of much research within a variety 
of fields, including optoelectronics and laser physics^, nonlinear optics|2], cavity 
quantum electrodynamics^, and biomolecular sensing 4 . Planar photonic crystal 
(PC) defect microcavities E], formed by an intentionally introduced perturbation 
to the periodic photonic lattice, have attracted attention for a variety of reasons, 
including lithographic control of many salient properties of the cavity modest E|, 
the potential for on-chip integration with PC waveguides 8 , and their ability to 
support resonant modes that simultaneously exhibit ultra-small (wavelength-scale) 
modal volumes and high quality factors^ E3j- In most of these works, a two- 
dimensional (2D) connected dielectric lattice of air holes is employed, often as a result 
of experimental constraints (for example, a connected lattice is necessary for undercut 
membrane geometries). The focus of these investigations has typically been on cavity 
modes with transverse electric (TE) polarization (magnetic field normal to the slab) 
in a hexagonal lattice, primarily because of the presence of a complete (in-plane) 
frequency band-gap for modes of such polarization^!]. Recent calculations^! E] 
and experiments JU|> however, have confirmed that the absence of a complete in-plane 
band-gap (in these cases, for TE modes within a connected square lattice of air holes) 
does not preclude the formation of low-loss, localized resonant states, provided that the 
modes are localized to regions of momentum space where the lattice provides strong 
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optical feedback. These results thus provide motivation for the study of localized 
transverse magnetic (TM) polarized (electric field normal to the slab) states within 
connected photonic lattices, which also lack a complete in-plane band-gap. 

Further encouragement is found in recent work incorporating a PC microcavity 
within a quantum cascade (QC) semiconductor heterostructure^3| to form an 
electrically-injected PC QC laser |15l 116) . The intersubband transitions responsible for 
photon generation within QC lasers are fundamentally TM in nature, while efficient 
electrical injection necessitates the use of a connected PC lattice. Though initial 
devices more employed more extended microcavity band- edge states, future efforts will 
involve optimization of current injection to create further localized electrically-driven 
defect mode lasers, attractive due to the potential for reduced cavity radiation loss 
(high quality factor) and further suppression of unwanted spontaneous emission. For 
such applications, an understanding of the salient properties of the TM defect states 
will be valuable. 

Our study involves two complementary techniques for understanding the 
properties of the TM-polarized defect states in square and hexagonal lattice 2D 
PCs. The first, detailed in Sectional is a symmetry-based analysis analogous to that 
considered in previous articles |Bl 1 17) for TE-like modes, where the tools of group theory 
are used to classify the symmetries and identify the dominant Fourier components of 
the defect states. 2D finite-difference time-domain (FDTD) simulations of the defect 
modes are presented in Section |3 and show a close correspondence with the group 
theory analysis, while complementing the results with a depiction of the real-space (in- 
plane) properties of the modes. The results of these sections provide a starting point 
from which more advanced defect designs can be generated, and a simple example of 
this is given in Section^ 

2. Symmetry Analysis 

The design of PC defect cavities has primarily been done through numerical methods 
such as FDTD. The overriding role that symmetry plays in physical systems suggests 
that the translational and rotational symmetries present in photonic lattices should 
play a leading part in determining the properties of the defect modes, and as a 
result approximate, less numerically intensive methods can also provide valuable 
information. In recent work|S| E|> the methods of group theory were used to 
analyze and classify the modes of 2D PC slab WG defect cavities for TE-like modes. 
The results of the symmetry-based analysis were compared with FDTD simulations 
and photoluminescence measurements, and showed a close correspondence between 
the approximate analysis, the detailed numerical simulations, and the experimental 
results. In this section, we extend the previously developed analysis to the case of 
TM-like modes, with a focus on QC laser applications. 

The discrete symmetries of the dielectric function, e(r), determine the symmetries 
present within Maxwell's equations. The group theory analysis we consider consists 
of two main steps. Modes of the unperturbed slab waveguide are used as a symmetry 
basis to generate approximate field patterns for the PC modes at the high symmetry 
points within the first Brillouin zone (IBZ) of the reciprocal lattice. The curvature 
of the photonic bands at these points are such that peaks and valleys are created in 
the energy-momentum dispersion surface. It is these peaks and valleys which are the 
seeds from which localized states are formed. The second step in our approach then 
utilizes these PC band-edge states as a basis to generate approximate forms for the 
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localized defect modes lying within the (partial) band-gap. 

We consider quasi-2D PC structures in this paper, where an effective index (n e ff) 
for the fundamental TM mode is used to take into account vertical waveguiding. In this 
quasi-2D system, the fundamental TM mode is represented by the scalar field E z ; in 
a true three dimensional (3D) PC slab WG system, modes can not be separated into 
pure TE and TM polarizations, but are nevertheless (for the lower lying frequency 
bands) significantly TE-like or TM-like in character. As a result, the approximate 
analysis presented here will have utility in studies of 3D systems as well (indeed, in 
previous work|fil IT7|. the close correspondence between the symmetry-based analysis 
and 3D FDTD results for TE-like modes was established). 

2.1. Hexagonal Lattice 



Figure 1. (a) Fundamental TM-like guided mode bandstructure (n c ff = 3.345) 
for a hexagonal photonic lattice with r/a=0.35. The air (cladding) light line 
is shown as a solid black line. The X and J-point regions from which donor-like 
defect modes are formed are circled, (b) Illustration of the (upper) real and (lower) 
reciprocal spaces of the two-dimensional hexagonal PC. The high-symmetry points 
of the hexagonal lattice, referenced to the center of an air hole, are a = (0,0), 
b = (a/2, 0), and c = (0, a/V3). 



The point group symmetry for TM modes in a 2D hexagonal lattice PC is 
Ce v ■ A plot of the approximate in- plane bandstructure (calculated by the plane- wave 
expansion technique^Hj) for the fundamental TM-like guided modes in this lattice is 
given in Figure ^a), where n c g = 3.345 in the dielectric regions, and the ratio of the 
hole radius r to lattice spacing a of r/a=0.35 has been chosen. 

For the fundamental TM modes of an unpatterned waveguide, the electric field 
within the mirror plane of the slab is given by Ek^ (rj_) = ze~^ ]l± ' r± \ where and 
are in-plane wavenumber and spatial coordinates, respectively. Upon patterning the 
slab waveguide, coupling occurs between waveguide modes with similar unperturbed 
frequencies and propagation constants modulo G, where G is a reciprocal lattice 
vector. This follows from the (approximate) conservation of frequency and crystal 
momentum. 

The high symmetry points of the first Brillouin zone (IBZ) of the 
hexagonal lattice (Figure ^| are the six degenerate A-points ({±(0, l)kx, 
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Table 1. Point Group character tables for the hexagonal lattice. 
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±(\/3/2, l/2)k x , ±(V3/2,-l/2)fcx}), the six degenerate J-points ({±(1/2, V3/2)kj, 
±(l/2,-V3/2)A:j, ±(l,0)/cj}), and the r-point=(0, 0), with fc x = and 
fcj = 47r/3a. The group of the wave vector, which defines the point group symmetry 
of a plane wave modulo G within the dielectric lattice, is C 2v , C^ Vl and Cq v for the 
X, J, and r points, respectively. Character tables^^l f° r these groups are given in 
Tabled 

Unlike the case for TE modes, the TM mode bandstructure does not exhibit a 
clear in-plane bandgap. Nevertheless, the curvature of the photonic bands is such 
that we can still expect localized defect states to form, provided that the modes are 
composed of Fourier components that lie in a region of momentum space where a 
bandgap does exist (i.e., a partial in-plane bandgap). In particular, the J-point of 
the third band (by frequency) has the correct curvature to form donor-type modes, 
and enough of a bandgap surrounding it to create localized states (the JT-point of the 
second band also has the correct curvature, but because the bandgap surrounding it 
is somewhat small we will not consider it). 

2.1.1. J-point Band-Edge Representation In determining an approximate description 
of the modes of the patterned slab waveguide we consider those modes of the 
unpatterned slab waveguide that are most strongly coupled considering temporal and 
spatial frequency conservation. At the J-point these modes have a scalar electric field 
(TM) given by E = ze~ lkj i' r± and a frequency uj j ps c\kj\/n e s. 

The star of k (*k) at the J-point, which consists of all the independent J-points 
within the IBZ, is given by (not uniquely) {k.j 1 , kj 2 }. A symmetry basis for the modes 
of the patterned slab waveguide at the Ji-satellite point can be found by applying the 
symmetry operations of the group of the wave vector (Gok., = C 3v ) to the seed vector 
Ek 7l - In this case, the basis is (Ek 7l , Ek /3 , Ek 7s ). Projecting this symmetry basis 
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onto the IRREP spaces of C 3v yields: 
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where Ai and £ label the IRREP spaces of C 3v (see Table Q] £ is a 2D IRREP), and 
the index a is used to denote the location of the origin at point a in the hexagonal 
lattice. Because the dielectric structure has a low index air hole at point a, we associate 

with the third highest frequency band mode. This is a result of the tendency, all 
things being equal, for modes with electric field concentrated within regions of high 
dielectric constant to be lower frequency than those with electric field concentrated in 
low dielectric regions [TT]. We associate the set of modes {E^i, E^}, with the first 
and second frequency bands, which, consistent with the bandstructurc calculations 
(Figure ^a)), are degenerate at the J-point. 

To fully define the modes at the J-point all modes of the *k must be included. 
The result of doing so is the following set of lower band modes, 
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and higher band modes, 
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The approximate J-point band-edge modes derived above have their origin at the 
center of an air hole (point a). The hexagonal lattice has two other high-symmetry 
points around which one may center a defect, points b and c shown in Figure ^b). 
Unlike point a, points b and c are of lower symmetry than that of the lattice (C2v and 
C 3Vt<7v> respectively). 

To determine the approximate band-edge modes for defect regions centered about 
these points, the lower and higher band modes are simply shifted to each of these 
points using the transformations = — b and rj_ = — c. The resulting bases 
are labeled J(l,2)f,, J3f,, J(l,2) c ,and J3 C for lower and higher bands centered about 
points b and c, respectively. 



2.1.2. J-point Donor Modes We now consider donor modes, which are formed by 
a local increase in the dielectric constant within the photonic lattice. As discussed 
previously, the lowest frequency region of k-space for which this can occur, while 
still enjoying the benefits of significant isolation from other bands, is the J-point 
of the third frequency band. Therefore, the appropriate symmetry basis to use for 
describing localized donor modes are the band-edge modes of J3 a , J3&, and J3 C for 
defects centered around points a, b, and c, respectively. For defects centered about 
point a the largest possible symmetry is that of the underlying hexagonal lattice, Cq v , 
whereas for defects about point b and point c the largest point group symmetries are 
C^v and C 3ViC r v , respectively. Correspondingly, the character values of representation 



K. Srinivasan and O. Painter, Design of two-dimensional photonic crystal. 



6 



Table 2. Symmetry and dominant Fourier components of the E-field for TM 



donor modes formed from the J3 band-edge in a hexagonal lattice. 
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ga, dj3 f J3 a Das i s under Cq v , representation S b,dj3 of the Jib basis under Ci v , and 
representation S 1 ^* 3 * of the J3 C basis under C^ Vi(Tv are given in Table From these 
character tables we find that these representations decompose as S a,dj3 = A'{ © B", 

gb, d J3 = Ai®Bi, and S c ' dj3 = E®A%. Using the appropriate projection operators [iH] 
on J3 a , a set of basis functions for the localized conduction band donor modes about 
point a is found. Similarly, projecting the basis functions J3f, onto the IRREPs of C 2v 
and the basis functions J3 C onto the IRREPs of C^ VtlJv , we get the donor modes about 
point b. The approximate forms of the donor modes found from these projections are 
given in Table [3 Also shown in this table is the new mode structure for a defect 
cavity of reduced symmetry (C 2v ) centered about point a (points b and c are already 
of lower point symmetry than the lattice), found by using the compatibility relations 
between the IRREPs of the full and reduced symmetry groups. 

The key pieces of information given by our approximate description are the 
dominant Fourier components that comprise each mode. Approximate forms for 
localized defect states can be produced by including some form of an envelope function 
along with these approximate symmetry-based forms; in previous work|20|. the effects 
of the underlying photonic lattice are captured through derivation of a Wannier-likc 
wave equation for the envelope of resonant TE modes. When combined with the 
results of the group theory based analysis, a very accurate picture (as confirmed 
by FDTD simulations) of the near-field behavior of the defect modes is produced. 
A similar Wannier-like equation for TM modes can be generated; however, for the 
purposes of this paper, we forego such an analysis and instead obtain more detailed 
information directly from FDTD simulations. The key point stressed in Ref. [201 that 
is of relevance here is that the ground state envelope functions transform effectively as 
the identity operator, and as a result have no bearing on the transformation properties 
(i.e., the classification of the symmetries of the modes) derived here. 

2.2. Square Lattice 

The point group symmetry of the square lattice photonic crystal is C^ v . A plot of the 
approximate in-plane bandstructure for the fundamental TM-like guided modes of a 
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square lattice of air holes with n e g = 3.345 for the dielectric regions and r/a=0.40 is 
given in Figure Ufa). 



Figure 2. (a) Fundamental TM-like guided mode bandstructure (n e g = 3.345) 
for a square photonic lattice with r/a=0.40. The air (cladding) light line is 
shown as a solid black line. The X and M-point regions from which defect modes 
are formed are circled in black, (b) Illustration of the (upper) real and (lower) 
reciprocal spaces of the two-dimensional PC with a square array of air holes. The 
high-symmetry points of the square lattice, referenced to the center of an air hole, 
are d = (0, 0), e = (0, a/2), and f = (a/2, a/2). 



The high symmetry points of the IBZ of the square lattice (Figure IHb)) are 
the four degenerate X-points ({±(1, 0)k x , ±(0, l)k x }), the four degenerate M-points 
((±v / 2/2,±\/2/2)fc M ), and the r-point=(0, 0), with k x = vr/a and k M = nV2/a. 
The group of the wave vector at the X-, M-, and T-points is Civ, C^ v , and C^ v , 
respectively. Character tables^!] f° r these groups are given in Table |3 

Based on the photonic bandstructure diagrams of Figure Ela), the candidate 
regions of k-space from which we expect to form defect modes include the X-point 
off the second frequency band (conduction band donor modes) and the M-point off 
the fourth frequency band (higher band donor modes). Discussion of the M-point of 
the first frequency band, also a candidate region of k-space (to create valence band 
acceptor modes) is not considered here in the interests of keeping the discussion as 
concise as possible. The analysis is quite similar to what was done for the hexagonal 
lattice; we begin by forming basis functions to describe the band-edge modes, and then 
couple these modes via the symmetry of the defect to get an approximate description 
of the cavity modes. 

2.2.1. X-point band-edge representation A symmetry basis for the modes of the 
square lattice PC at the X-point is found by applying the symmetry operations of 
the group of the wave vector (£7 k x = Civ) to the seed vector Ek Xl . In this case, the 
basis is simply (Ek Xl , E_k Xl )• Projecting this symmetry basis onto the IRREP spaces 
of Civ , and then including all the modes of the star of *k, results in the following set 
of degenerate valence band modes, 
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Table 3. Point Group character tables for the square lattice. 
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and degenerate conduction band modes, 

X2 d = I ( C ° S ^ ' % , (5) 
Vcos(k X2 • r±)J 

In the square lattice there are three different high-symmetry points around 
which one may center a defect, labeled d, e, and / in Figure EJb). Points d and 
/ maintain the C^ v symmetry of the lattice, while point e has the lowered symmetry 
C'iv,<j v ■ The shifted valence and conduction band-edge modes are found through the 
transformations r"^ = rx — e and = — f, and are labeled as Xl e , X2 e , XI f, 
and X2f for points e and / respectively. 

2.2.2. M-point band- edge representation The group of the wavevector at the M-point 
is C± v , and as such the symmetry basis, S M , includes all the M-points of the IBZ, 
S M — (Ek M ,Ek M2 ,E_k Ml ,E_k M2 ). As determined from its character under C± v 
(Table©, S M = E®A'{®B 2 \ The doubly degenerate IRREP E represents the second 
and third frequency bands as they are degenerate at the M-point as shown in Figure 
HJa), while the B% and A'{ IRREPs represent the first and fourth frequency bands, 
respectively, as a result of the lattice being centered on an air hole. Approximate 
forms for the modes of these bands at the M-point are calculated by projecting the 
symmetry basis onto the three IRREP spaces. With the origin centered at point d, 
the band-edge modes are: 

M l d = z (cos(k Ml • ri) - cos(k M2 ■ r^)) , (6) 

MA d = z (cos(k Ml • rl) + cos(k M2 ■ . (8) 

As was done for the X point, shifted band-edge modes for points e and / are easily 
found and are labeled Ml e , M(2,3) e , M4 e , and Ml/, M(2,3)/, and M4/. 
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Table 4. Symmetry and dominant Fourier components for the E-field of TM 
donor modes formed from the X2 and MA band-edges in a square lattice. 
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2.2.3. Donor modes For the square PC lattice, as described earlier, two band- 
edge minima from which donor modes can be formed are the X-point in the second 
frequency band and the M-point in the fourth frequency band (Figure Ufa)). For 
a symmetric defect centered at points d and / the band-edge modes at the X 
and M points couple to form resonant modes which transform as IRREPs of C± v , 
whereas for the lower symmetry point e the defect modes transform as IRREPs 
of C 2 v,a v - The representations describing how the X2d, X2 e , and X3f symmetry 
bases transform under the appropriate point group are given by S d ' dx2 , S e ' dx2 , and 
S*' X2 , respectively. Similarly, the representations describing how the MAd, MA e , 
and MAf symmetry bases transform under the appropriate point group are given by 
S d ' dM4 , S e ' dMA , and S f ' dM4 , respectively. From their characters in Table we find 
that S d - dx2 = A'l 8 B", S e4x2 = ii ffi B 2 , and S f ' dx2 = E, while S d - dM4 = A", 
ge,d M 4 — anc j gf,d M 4. _ p or p i n t s d and /, defects may also be formed with 

lower symmetry than the Ci v symmetry of the square lattice. Using the compatibility 
relations between the IRREPs of the full and reduced symmetry groups, we determine 
the new mode structure for a defect of C% v cr„ symmetry (i.e., with mirror planes 
along the x and y directions of Figure |2Jb)). Projecting the different symmetry bases 
onto the IRREPS of the point group symmetry of the different cavities one obtains 
approximate forms for the TM donor modes in the square lattice, the results of which 
are tabulated in Tabled 

|| The symmetries and fundamental momentum components of the possible donor modes formed from 
the MA band-edge are trivially given by a single band-edge mode. This is due to the fact that the 
Af-point in the square lattice is highly symmetric, and the group of the wave vector mixes all four of 
the M-points of the IBZ. 
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3. FDTD Simulations 

In this section, FDTD simulations of the donor defect modes identified in Section|21arc 
presented. The FDTD calculations establish the effectiveness of the chosen symmetry 
bases of the group theoretical analysis in describing the defect cavity modes. To reduce 
computational time, we limit our study to 2D FDTD simulations with an effective 
index used to account for vertical guiding; previous papers jBJ El have shown that 3D 
FDTD simulations of TE-like modes in full PC slab WG defect cavities correspond well 
with the modes found through the approximate group theory analysis. The FDTD 
calculations were performed on a mesh with 20 points per lattice spacing. Excitation 
of the cavity modes was performed by an initial field (E z ) with a localized Gaussian 
profile, located in a position of low symmetry so as not to exclude any possible resonant 
modes. TM modes are selected out by constraining the normal component of the field 
to be purely E z , and a pair of mirror planes (cr x , o~ v ) were used to filter out cavity 
modes according to their projection on the IRREPs of Ci v . Each FDTD-generated 
cavity mode is thus labelled by the C2 V IRREP by which it transforms and an index 
corresponding to its energy (frequency) level. 

FDTD-generated results are compared with the results of the symmetry analysis 
both through confirmation of the symmetry classifications under the IRREPs of Ci v 
(although a careful analysis of all available mirror planes within the FDTD simulations 
would allow classification under higher symmetries such as C± v or Cq v ) and through 
an analysis of the dominant Fourier components of the modes (through examination 
of the spatial Fourier transform of E z (E z )). 

3.1. Hexagonal Lattice 



We begin by considering donor-type defect modes in the hexagonal lattice formed 
around the J-point of the third frequency band, where from our symmetry analysis 
we expect to see modes of A\ and B\ symmetries (under C% v ). In Figure we show 
two examples of simple defects within the hexagonal lattice, one localized around the 
point a of the lattice, the other around point b. The chosen defect around the b- 
point in the hexagonal lattice, the defect cavity we will consider here, consists of two 
reduced size air holes along the y-axis, with r/a = 0.35 and r<j f/a = 0.30 (Figure 
EJb)). The spectrum for A\ modes shown in Figure a) shows a single dominant 
peak in the frequency range of interest; the real space and Fourier space versions 




r/a 



r/a=0.35 



Figure 3. Defect geometries for donor modes about the (a) a-point and (b) 
b-point. 
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of this field are given in Figure Elb)-(c), and indicate that the mode has dominant 
Fourier components at the six J-points within the IBZ, as predicted by our symmetry 
analysis. The Fourier transformed electric field for the B\ symmetry mode is shown in 
Figure^f); the mode has its dominant Fourier components at the six J-points, and is 
somewhat delocalized in comparison to the A\ mode (FigureQJe)). Similar results and 
correspondence were found for the defect around point a of the lattice (Figure Of a)). 
As a final comment, we note that FDTD-generated spectra for A 2 and B2 symmetry 
modes were also examined (for both points a and b), and significant resonance peaks 
were not seen in the frequency range of interest. 
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Figure 4. Donor modes about the b-pt. in the hexagonal lattice, (a) Spectrum 
for A\ symmetry modes, (b) E z and (c) \E Z \ for A\ mode, (d) Spectrum for B\ 
symmetry modes, (c) E z and (d) \E Z \ for B\ mode. 



3.2. Square Lattice 

In Section|2]wc identified three high symmetry points in the square lattice about which 
defect modes are formed. For each of these points, we considered donor modes off the 
X point of the second frequency band and M point of the fourth frequency band. We 
now present FDTD simulations which verify the symmetry classification and Fourier 
space composition for the defect modes predicted by the group theory analysis. 

Figure [S] shows three examples of defect cavities formed about the d, e, and 
/ points of the square lattice. We present here the FDTD results for the defect 
cavity about the / point. In this cavity a localized defect state is formed that has z- 
component of the electric field predominantly within the high refractive index material, 
which for application as a QC laser cavity would provide the largest possible overlap 
factor with the gain region in the semiconductor. Donor modes about the /-point 
are created by reducing the size of four air holes (Figure E{c)). For this example, 
we have taken r/a — 0.40 and rdcf/a = 0.35. The symmetry classifications for the 
modes predicted in SectionElare confirmed as this perturbation to the dielectric lattice 
creates donor modes of type B\, B2, and A2 symmetry, with the primary B\ mode 
at a normalized frequency of a/Ao = 0.220 (Figure EJa)), the dominant B2 mode at 
a normalized frequency of o/Xq — 0.220 (Figure Eld)), and the main A2 mode at 
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Figure 5. Defect geometries for donor modes about the (a) d-point and (b) 
e-point and (c) /-point. 
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Figure 6. Properties of /-pt. modes in the square lattice, (a) B± symmetry mode 
spectrum, (b) E z and (c) \E Z \ for low frequency B± mode, (d) B2 symmetry mode 
spectrum, (e) E z and (f) \E Z \ for low frequency B2 mode, (g) A2 symmetry mode 
spectrum, (h) E z and (i) \E Z \ for high frequency A2 mode. 



a/Xo = 0.419 (Figure Etg))- The real and Fourier space versions of the electric field 
for the low frequency peaks show that the B\ mode (Figure EJb)-(c)) has dominant 
Fourier components at ±X2, and the B2 mode (Figure Ete)-(f)) has its dominant 
Fourier components at ±Xi, consistent with the group theory analysis. These modes 
are in fact degenerate (consistent with the group theory prediction that they originate 
from a 2D IRREP for a defect of C& v symmetry) and act as a dipole-like pair of 
modes, in analogy with the hexagonal lattice dipole modes studied in previous work 6 . 
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Finally, |E| and |E Z | for the higher frequency A 2 peak are shown in Figure EHh)-(i), 
and show that it is composed of dominant Fourier components at the four M points. 

4. Summary and Discussion 

The group theory analysis presented in Section can be used to ascertain a number 
of important properties about TM-polarized defect modes in connected PC lattices, 
and as such, serves as a starting point from which more careful designs of PC defect 
cavities can originate. For example, recognizing potential modal degeneracies and 
modifications that can be made to split or remove these degeneracies can be important 
in laser design[7| 0. Coupling of the cavities to other in-plane elements can benefit 
from an understanding of the modes' dominant in- plane emission directions, while 
knowledge of the position of field nodes/antinodes with respect to the dielectric and air 
regions is important as some applications (such as coupling to semiconductor quantum 
dots) benefit from field maxima in the dielectric regions, while others (such as coupling 
to states of cold atoms "dropped" within the cavity) are better served by field maxima 
in the air holes. 

With the results obtained from the group theory analysis in hand, more 
sophisticated modeling, such as full 3D FDTD simulations, can be used to determine 
more quantitative properties of cavity modes such as modal frequency, cavity quality 
factor (Q), and precise near and far-field emission pattern. This approach was used in 
Ref. which focused on the design of high-Q wavelength-scale PC microcavities. In 
this work, the group theory based approach was used to analyze TE-polarized modes 
localized to different high symmetry points within connected square and hexagonal 
lattice PCs. Fourier space design principles were then used to determine which of 
these candidate modes were particularly suitable for high-Q applications, and then 
full 3D FDTD modeling was used in the final design steps. 




Figure 7. Properties of Jf-point E^ x2 symmetry mode in a graded square 
lattice, (a) E z (b) \E X \. 



As a simple example, we conclude by considering the formation of a PC defect 
laser cavity for TM-polarized modes (used in conjunction with QC heterostructures, 
for example) . The square lattice defect cavities considered in Section [21 had fairly 
large-size air holes (r/a — 0.45), which in fabricated devices can be problematic for 
a number of reasons, including reduced dielectric-field overlap (r) and poor thermal 
stability due to the porosity of the lattice. In electrically injected devices, this porosity 
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could also hinder defect mode lasing, as pumping of the central (defect) region of the 
cavity will be more difficult for such small current paths. On the other hand, larger 
holes create larger frequency-space gaps in which the modes can be localized. One 
method for achieving high in-plane reflectivity while reducing the porosity of the lattice 
is to employ a graded lattice design^! ■ Starting with the requirement that our lasing 
mode have an electric field antinode in a dielectric region (this mode is similar to the 
E^^ 2 mode of the /-point defect cavity), we use our group theory based analysis 

to select the low- frequency ¥fj^ X2 mode as a mode of interest. We then employ the 
graded lattice geometry shown in Figure Ufa) to create such a mode. It consists 
of two levels of confinement: first, an initial jump in hole radius between the holes 
centered at (0, ±o/2) (which have r/a = 0.225) and their nearest neighbors (which 
have r/a — 0.25), and then, a parabolic increase in hole radius moving outwards from 
the defect center (with a maximum value r/a — 0.30, equal to the hole size used in 
lasing devices in Ref. dO). The initial jump in r/a serves as a potential well to 
help confine the mode in real space, while the grade allows for the hole radius to 
be kept small in the region where the mode is primarily located (to keep T high), 
but large enough outside of this region to produce a high in-plane reflectivity. The 
resulting real space field E z and Fourier transform \E Z \ for the Ai mode are shown 
in Figured Based on these images, we see that our initial objectives have been met: 
the mode remains well localized to the central region of the cavity, despite having now 
chosen considerably smaller holes. With this lattice geometry taken as a template, the 
Fourier space design principles of Ref. ^] could now be applied, in conjunction with 
3D FDTD simulations, to optimize the cavity Q factor for use in high performance 
devices. 
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